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We show that in the case of the BFKL pomeron with running coupling 
the diffusion pattern is strongly modified and is characterised by the sudden 
\Q , tunneling transition to the non-perturbative regime. We suggest that by 

using the b - expansion method one can suppress the non-perturbative 
Pomeron and isolate purely perturbative part of the gluon Green's function. 
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PACS numbers: 12.38 Cy 

One of the major issues in the high energy limit of QCD is the problem 
of interplay between the perturbative and non-perturbative regimes. The 
BFKL Pomeron in the leading logarithmic approximation [1] results in the 
equation which exhibits the characteristic diffusion pattern [2]. In the pro- 
cess with two, comparable hard scales (for example 7*— 7*, or forward jet /tt° 
in DIS) the distribution of the transverse scales of the gluons broadens with 
increasing rapidity, with the width At ~ y/c^Y where t — In A: 2 /A 2 with k 2 
rS ' being the transverse momentum of the gluon and A being a QCD parameter. 

Thus, for sufficiently high energies the distribution of the gluon momenta 
will always reach non-perturbative regime. This picture of diffusion is well 
rooted in the case of the fixed a s coupling. The effect of subleading cor- 
rections [3] is, among other, the running of the QCD coupling. In this case 
one would expect small modification of the diffusion picture, namely that 
the distribution of the momenta will develop an asymmetry towards lower 
scales. The exponent of the cross section gains additional term b 2 a^Y 3 apart 
from the usual leading a s Y term, where Y is the rapidity of the process. 
This picture holds however for rapidities which are not too large. In [4] it 
has been pointed out that the transition to the non-perturbative regime can 
occur as a sudden tunneling effect rather than a gradual diffusion process. 
Instead of slow increase of the width of the distribution of the momenta, one 
observes a sudden transition from the perturbative scales t r ~ t = lnQ 2 /A 2 
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into the non-perturbative regime t' ~ t > 0, where i — lnQo/A 2 is the scale 
set by the regularisation of the running coupling. The characteristic feature 
of this effect is that it occurs without any passage through the intermediate 
scales 1 . 

To gain insight into the effect of the tunneling we shall consider the 
small x evolution which is controlled by the BFKL equation of the form 

?^-=K 9 G (!) 

where K is the BFKL kernel and G(Y, t,to) gluon Green's function eval- 
uated at rapidity Y = lnl/x and scales t and to with initial condition 
G(Y = 0;t,to) = S(t — to). Kernel K = a s K Q is the usual BFKL kernel in 
the leading logarithmic approximation, but we additionally introduce the 
running of the coupling a s (t) . 

The detailed mechanism of the tunneling is illustrated in Fig.l where we 
show the contour plots in the t and Y plane of the function 

!(Yy . t () _ gm^-^m (2) 

which illustrates the change of distribution of the transverse momenta dur- 
ing evolution between two points: from (0, to) to (Y,t). At the beginning 
of the evolution Y < 50 the solution exhibits typical diffusion pattern 
and the broadening of the distribution. Around Y = 60, 70 the tunnel- 
ing transition occurs in which emerges second region concentrated around 
non-perturbative scale t. At higher rapidities the evolution is concentrated 
only in the non-perturbative region. One can estimate the value of rapidity 
at which the tunneling transition takes place. We just have to compare the 
contribution to the Green's function which comes from the tunneling con- 
figuration with that coming from the perturbative evolution. The tunneling 
configuration can be written as: 

G tunnel (y; t, t) ~ e-M^V^/ 2 (3) 

where the first term corresponds to branching from t to t, second to the 
evolution at t with a non-perturbative Pomeron intercept ujp and the last 
term corresponds to branching from t back to t. Instead the configuration 
which comes from the normal perturbative evolution is just 

Gp_ert(^; t, t) ~ e^ Y (4) 

1 It has to be stressed that the unitarity effects can in principle change significantly 
the phenomenon of tunneling. It was noticed in [5] that in the case of the non-linear 
small x evolution equation, the generation of the saturation scale Q s (x) leads to the 
suppression of diffusion into the low scales k < Q s (x) and the distribution of the 
momenta is driven towards the perturbative regime. 
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where uj s (t) is the effective perturbative saddle point exponent. 
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Fig. 1. Contour plots for f(Y,y,t,t f ), illustrating different stages of the evolution: 
Y = 30 and 50 illustrates standard 'cigar' type plots, Y = 60-70 show the point 
where tunneling begins to play a role, while for Y > 80 the momentum configura- 
tion is only in the non-pert urbative regime. 



Tunneling transition will take place when Gtunnei — G per t and this occurs 
at rapidity 

^tunnel (t) = y-r (5) 

which shows the linear dependence on t and the slope governed by the value 
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of the intercept of the non-perturbative Pomeron cjjp. 

From this analysis it is clear that in the case of the small x evolution with 
running coupling one has always contamination from the non-perturbative 
contribution. The real problem lies then in the definition of the perturbative 
hard Pomeron. A way of solving this problem is to consider the 6-expansion 
[6] i.e. take the limit in which b ^ where b is the beta function coefficient 
of the QCD coupling. It has been shown [6] that the gluon Green's function 
can be decomposed into the perturbative and non-perturbative terms and 
that the ratio of these two components is roughly of the form 

where g(a s (t)) is some function which can depend on the details of the model 
for small x evolution (for example Airy, collinear, full BFKL). We see, that 
the non-perturbative Pomeron is asymptotically leading since we always 
have ujjp > uj s (t) but is suppressed by the universal exponential factor. By 
taking the limit b — > we can eliminate the non-perturbative Pomeron and 
are able to isolate purely perturbative contribution which is independent of 
the given regularisation procedure for the running coupling. The b ^ limit 
corresponds to the assumption of the very slowly varying QCD coupling 
(expansion around the fixed coupling limit). Using the b - expansion, see [6], 
one is able to identify systematically various types of diffusion corrections 
~ 6 2 a^y 3 , b 2 a^Y 2 and specify the radius of convergence for this series given 
by the parameter ( c = bxm&s 2 (to) — 0.264, where Xm = 4 In 2 being the 
minimum of the BFKL kernel eigenvalue. We have also checked numerically 
that the perturbative behaviour indeed breaks down at rapidities Y ~ t 2 
which can be seen from the above quoted form of the parameter £ c . In 
Fig. 2 we show the maximum rapidity for which the perturbative part can 
be defined using two methods. In the first one, we calculate the solution 
to the equation (1) with two different regularisations for a 8 (t) and define 
y max as the limiting rapidity at which the two solutions start to diverge - 
dashed line. In the second method we apply the 6-expansion method, and 
calculate the perturbative gluon Green's function from the following series 
lnG P ert(^; t) = J2i b l li(a s (t),Y). We then truncate this series at different 
orders and define again Y max as the divergence point. We clearly see from 
Fig. 2 that the first method yields linear dependence on t which is consistent 
with formula (5) and suggest the tunneling as the relevant mechanism for 
the breakdown of perturbative evaluation. In the second case, we observe 
~ t 2 behaviour which is consistent with the prediction based on the critical 
value ( c . Thus the perturbative prediction can be in principle made up to 
the values Ym ax ~ t 2 . 
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Fig. 2. The maximum perturbatively accessible value of Y, as a function of l/a 8 (to), 
determined by comparing different non-perturbative regularisations of a s , or dif- 
ferent truncations of the 6-expansion. 
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